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ABSTRACT 

We construct a new model of Type la Supernovae (SNe la), based on the 
single degenerate scenario, taking account of the metaUicity dependences of the 
white dwarf (WD) wind and the mass-stripping effect on the binary companion 
star. Our model naturally predicts that the SN la lifetime distribution spans a 
range of 0.1 — 20 Gyr with the double peaks at ~ 0.1 and 1 Gyr. While the 
present SN la rate in elliptical galaxies can be reproduced with the old popula- 
tion of the red-giants+WD systems, the large SN la rate in radio galaxies could 
be explained with the young population of the main-sequence-|-WD systems. Be- 
cause of the metallicity effect, i.e., because of the lack of winds from WDs in the 
binary systems, the SN la rate in the systems with [Fe/H] < — 1, e.g., high-z 
spiral galaxies, is supposed to be very small. Our SN la model can give better 
reproduction of the [{a, Mn, Zn)/Fe]-[Fe/H] relations in the solar neighborhood 
than other models such as the double-degenerate scenario. The metallicity effect 
is more strongly required in the presence of the young population of SNe la. We 
also succeed in reproducing the galactic supernova rates with their dependence 
on the morphological type of galaxies, and the cosmic SN la rate history with 
a peak at 2; ~ 1. At z > 1, the predicted SN la rate decreases toward higher 
redshifts and SNe la will be observed only in the systems that have evolved with 
a short timescale of chemical enrichment. This suggests that the evolution effect 
in the supernova cosmology can be small. 

Subject headings: galaxies: abundances — galaxies: evolution — stars: super- 
novae 
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1. Introduction 



The star formation liistories of galaxies are imprinted on stellar populations of the 
galaxies. Because of the age-metallicity degeneracy for stellar populations, it is difficult to 
determine "age" from their colors and spectra alone. However, the elemental abundance 
ratios can provide independent information on "age", because differ ent heavy e l ements are 
prod uced from different supernovae with different timescales (e.g., iPagell 119971 : iMatteucci 
20011 ). To derive the information of "ages" from elemental abundance patterns as a cosmic 
clock, the most important uncertainty is the lifetime of Type la Supernovae (SNe la), i.e., 
the evolutionary time (delay-time) of the progenitor from the main-sequence through the 
explosion, which is mainly determined by the lifetime of the relatively low-mass (< 8Mq) 
companion star of the exploding white dwarf (WD). 



There exist two distinct types of supernova explosions (e.g., lArnettI Il996l : iFilippenko 



19971 ): One is Type II supernovae (SNe II), which are the core collapse-induced explosions of 
massive stars (> 8Mq) with short lifetimes of 10^~^ yrs, and produce more a-elements (O, 
Mg, Si, S, Ca, and Ti) relative to Fe with respect to the solar ratios (i.e., [a/Fe] > 0). Type 
lb and Ic Supernovae (SNe lb and Ic) are also the core-collapse supernovae, but have lost 
the envelope material before the explosions. Recently, it is found that brigh t core-collapse 



super novae, hypernovae (HNe), produce an important amount of Fe (e.g., iNomoto et al. 
20061 ). In chemical evolution models, the contributions of SNe lb, Ic, and HNe can be 



included in that of SNe II (iKobayashi et al.ll2006l ). The other is SNe la, which are the 
thermonuclear explosi ons of accreting WDs in close binaries and produ ce mostly Fe and 



little a-elements (e.g., iNomoto et al.lll994l : iHillebrandt &: Niemeyerl l2000l ) . The progenitors 
of the majority of SNe la are most likely the Chandrasekhar (Ch) mass WDs (e.g., Nomoto, 
Iwamoto & Kishimoto 1997 for a review; see also Nomoto et al. 2007), although the sub-Ch 
mass models might correspond to some peculiar subluminous SNe la. The early time spectra 
of the majority of SNe la are in excellent agreement with the synthetic spectra of the Ch 
mass models, whil e the spectra of the sub-Ch mass models are too blue to be compatible 
with observations (IHoflich fc Khokhlovl 1 19961 : iNugent et al.l 119971 ) . 



For the evolution of accreting C+0 WDs toward the Ch mass, two scenarios have been 
proposed; One is a double-degenerate (DD) sce nario, i.e., merging of double C-l-0 WD s with 
a combined mass surpassing the Ch mass limit (llben fc Tutukovlll984l : IWebbinklll984l ). From 
the long-term observational search (SPY project) of the WD binaries for the DD systems, 
more than thousand WDs are found. However, their combined masses are smaller than the 



( Napiwotzki 


2007; 


Geier et al. 


2007 



Tovmassian et al.ll2008l ). Theoretically, it has been suggested that the me r ging o f WDs leads 



to accretion-induced collapse rather than SNe la (ISaio fc Nomotdll985l . Il998l ). The other 
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is a single-degenerate (SD) scenario, i.e., the WD mass gro ys by accretion of hydrogen-rich 



matter via mass transfer from a binary companion (e.g., iNomoto et aLll2000l : iLivid 12001 
for reviews). Two progenitor sys tems have been found for red-giant (RG) and near main- 



sequence (MS) companions (e.g., iHachisu et al.lll999bl ). The SD scenario has predicted the 



detection of hydrogen-rich circumstellar matter from the companion star for a certain class 



Hamuv et al. 


2003; 


Pat at et al. 


2007) 



The lifetime of SNe la has been estimated from different aspects as follows, and these 
estimates conflict each other; (1) Population synthesis — The lifetime distribution of SNe la 
has b een predicted with the population synthesis models for the evoluti on of binary systems 
(e.g. , lYungelsonll2005l ; iBelczynski et al.ll2005l ; iTutukov &: Fedorovall2007l ) . The predicted rate 



of SD systems in those models is much smaller tha n DD systems, and the majority of SNe 
la have the lifetime as short as ~ 0.1 Gyr (e.g., iRuiz-Lapuente. Burkert. fc Canall Il995 



Yungelson fc Livid Il998l ). However, the optically thick winds from the WD, which were 
not taken into account in the above models, have been shown to play an essential role in 



the ev olution of the accreting WDs (jHachisu. Kato fc Nomotd Il996l . Il999al ; iHachisu et al. 



1999bl ). Furthermore, the strippirig effe ct may increase the parameter space for the SD 
systems (IHachisu. Kato fc Nomotd l2008l ). If these effects are taken into account in the 
population synthesis models, the DD rate could be smaller, and the SD rate could be larger. 

(2) Chemical evolution — The typical lifetime of SNe la could have been constrained 
from the chemical evolution of the Milky Way Galaxy, most notably in the [a/Fe]-[Fe/H] 
relations. Metal- poor stars with [Fe/H] < — 1 have JQ/Fej ~ 0.5 for Mg, Si, and Ca on 
the average (e.g.. iMcWilliam et al.lll995l ; ICayrel et al.ll2004j ). whil e disk stars with [Fe/Hl 
> — 1 show a dec rease in [a/Fe] with increasing metallicity (e.g., lEdvardsson et al. 



1993 



Bensby et al.ll2004j ). Such an evolutionary change in [a/Fe] against [Fe/H] has been explained 
wit h the early heavy-elernent p roduction by SNe II and the delayed enrichment of Fe by SNe 
la (iMatteucci fc Greggi ol ll986h . Conversely, chemical evolution models can constrain the 
nature of progenitor systems of SNe la; for example, Yoshii, Tsujimoto, & Nomoto (1996) 
estimated the lifetime of SN la progenitors to be as long as 0.5 — 3 Gyr. 

(3) Supernova rates — The redshift evolution of the SN la rate can give a rough estimate 
of the typical lifetime of SNe la (i.e., delay-time) for the assumed star formation history. 
The decrease of the SN la rate seen at z ~ 1 corresponds to the delay-time of ~ 3 Gyr 



(jStrolger et al.l 120041 ). However, the observed cosmic star formation rates (SFRs) involve 
uncertainties from dust extinction and completeness. The environmental dependence on 
galaxy evolution must be important if the SNe la rate depends on metallicity. SN la rates 
depending on the galaxy type can give stronger constraint on the SN la progenitor models 
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( ICalura fc Matteuccil l2006l ). iDella Valle et al.l (120051 ) suggested that the short hfetime of 
~ 0.1 Gyr is r equire d to explain the high SN la rate in radio-loud elliptical galaxies, and 
Sullivan et al.l (120061 ) suggested the similar lifetime distribution from the relation between the 
SN la rate and the specific SFR. Therefore, th e bimodal distribution of the lifetime of SNe la. 



so-ca lled "A-l-B" model, has been proposed (IScannapieco fc BildstenI l2005l : iMannucci et al. 
2006h . 



In addition, Kobayashi et al. (1998, hereafter K98) have found the metallicity effect on 
the WD winds, consequently, on the occurrence of SNe la, based on the simulations of the 
progenitor evolution by Hachisu et al. (1999ab). If the iron abundance of the progenitors is 
as low as [Fe/H] < — 1.1, then the wind is too weak for SNe la to occur. Thus the SN la 
rate in the system with [Fe/H] < — 1.1 is supposed to be "very small". This prediction has 
not been observatio nally confirmed yet, and several S Ne la in metal-poor environment have 
been reported (e.g., iPrieto. Stanek. fc Beacomll2008l ). Statistical studies of a large sample 
are required, since there must be a spread in the metallicity distribution function of stars in 
a galaxy and also the [0/Fe] is uncertain (see §2.31 for the details). 

If the lifetime distribution function of SNe la depends on metallicity, the detailed cal- 
culation of chemical evolution is even more important to understand the observations of 
supernova rates and abundance ratios. The redshift evolution of the cosmic SN la rate 
is affected by chemical enrichment histories in various types of galaxies. The estimates 
of the star formation histories of galaxies from their abundance ratios are affected by the 
metalhcity effect on th e SN la rate. In addition, Fe can be produced by the four sources 
( iKobayashi et al.l 120061 ): i) SNe la, ii) relatively large Fe production from low-mass SNe II 
with M = 13 — 15Mq, iii) large Fe production from HNe, or iv) the SN 1.5 explosion of 
some AGB stars. Therefore, not only [a/Fe] but also other elements should be considered 
to discuss the contribution of SNe la. 

This may be the case for dwarf spheroidal (dSph) galaxies; the observed low [a/Fe] 
at low [Fe/H] in dSphs is often su ggested to be du e to the significant contribution of 



SNe la at such low metallicity (e.g., I Venn et al.l l2004j ) . However, there is a large variety 



among galaxies, and a decreasin g trend in [M g/Fe] at [Fe/H] < — 1 may be seen in Draco 
( IShetrone. Cote fc Sargentll200ll ) and Sculper (IShetrone et al.ll2003l ). but not for other a ele- 



ment s such as Si and Ca . For more metal-poor stars, so me stars show high [g/Fe] (iKoch et al. 



20081 : iFrebel et al.ll2009l ). but others show low [a/Fe] (lAoki et al.ll2009l : iFrebel et al.ll2009tf ). 
Since [Mn/Fe] in dSph galaxies is as low as in Milky Way ha lo stars, the abundance pattern 
can be better explained with low-mass SNe II than SNe la (ITolstoy et al.ll2003l : IKobayashi 



2003 



Kobayashi et al.l l2006l . hereafter K06). This is also supported by the low [Ba/Fe] in 



Aoki et al.l (120091 ). Chemodynamical evolution of dSph galaxies are complicated because 
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of the low SFR and large dark matter content, and can be largely affected only by a few 
supernovae. At least, the low [a/Fe] in low metallicity stars does not necessarily mean the 
contribution of SNe la, and certainly more detailed studies are required under the reasonable 
model of SN la progenitors. 

Similarl y, it is argued that damped Lyman a (DLA) systems show low [a/Fe] at low 



[Fe/H] (e.g., IPettini et al.lll999[ ). IProchaska fc Wolfd (l2002f ). how ever, concluded that [Si/Fe] 



0.3 at [Fe/H] < —1.5, and also similar values were presented in lDessauges-Zavadsky et al. 



(120071 ). The elemental abundance ratios in DLA systems are affected by the dust depletion. 
The dust depletion has been often estimated using [Zn/Fe] ratio, but there is no assur- 
ance that [Zn/Fe] ~ for low-metall icity. In fact, an in c reasing trend of [Z n/ Fe] has been 
observed in the Milky Way Galaxy JPrimas et al.l I2OO0I : iNissen et al.l 120071 ). iNissen et al.l 



(I2OO7I ) showed that [S/Zn], which is free from the dust effect, of halo stars is similar to that 
in DLA systems. Such [S/Zn] and [Si/Fe] could be explained by hypernovae. Since the 
observation is for neutral gas in galaxies, it can be eas ily affected by dilution of hydrogen 
in chemo dynamical simulations (IKobayashi et al.l 120071 ) . Again, more detailed studies are 
required. 

The modeling of SNe la in the chemical evolution was first made by Greggio & Renzini 
(1983), where the SN la rate is calculated from the distribution of binaries based on Whelan 
& Iben (1973), and then extended by Greggio (1996) including the efficiency of acc retion of 
the en velop e of the secondary to t he WD. This model has recently been updated in iGreggio 
(I2OO5I ) and iMatteucci et al.l (120061 ). An alternative model was proposed by K98, fully based 
on the SD scenario, where the metallicity effect is essential to explain the [0/Fe]-[Fe/H] 
relation in the solar neighborhood. Kobayashi, Tsujimoto, & Nomoto (2000, hereafter KOO) 
applied this model for different types of galaxies, and succeeded in reproducing the present 
SN II and la rates in different types of galaxies. Matteucci & Recchi (2001, hereafter MROl), 
however, argued that they could not reproduce K98 results and that the lack of SNe la at low 
metallicities produces results at variance with the observed [0/Fe]-[Fe/H] relation. However, 
this is due to their misunderstanding of the parameterization of the K98 models, as described 
in detail in §5. 

In this paper, we first show the lifetime distribution functions derived from different SN 
la models. We then introduce ou r new SN la model by taking into account the stripping effect 
(IHachisu. Kato fc Nomotdl2008l ) and show the lifetime distribution function as a function of 
metallicity (§2). We then evaluate the SN la models from the comparison with observations. 
The evolution of the elemental abundance ratios in the solar neighborhood (§3) and the 
galactic and cosmic supernova rates (§4) can be better reproduced with our SN la model 
including the metallicity effect. In §5, we discuss the two formulations in K98 and MROl in 
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detail, and reproduce both results using the alternative formulation. In §6, we summarize 
our conclusions. The detailed discussion on the elemental abundance ratios in other galaxies 
will be in the second paper. 



2. Lifetime Distribution 



2.1. Previous Models of Type la Supernovae 

Lifetime (i.e., delay-time) of SNe la is not a single parameter. The SN la lifetime, 
tia, can be determined from the lifetime of companion stars in C+O WD binary systems, 
because the time interval from the start of mass accretion to the explosion is as short as 
< 10^ yr. Since the companion's mass ranges are different among various SN la scenarios, 
the typical lifetime is different for the different SN la models. In chemical evolution models, 
two formulations have been proposed in KOO and MROl (see §5 for details, Eqs.[2] and [1]). 
For the DD scenario, the timescale of gravitational radiation until the two WDs start to 
merge is also important, and the lifetime distribution functions have been provided from the 
population synthesis of the binary evolution. 

Figure [1] shows the distribution functions of the companion's mass ma (upper panel) 
and lifetime tia (lower panel) for the three SN la models. These correspond to the SN la rate 
of the simple stellar population, which is defined as a single generation of stars with the same 
age and metallicity. We adopt the Salpeter initial mass function (IMF) (j){m) oc m~^'^^ for 
m = 0.07 — SOMf;^, and the met allicity dependent main-sequence lifetimes, w hich is calculated 



bv iKodama fc Arimotd (119971 ) with the stellar evolution code described in llwamoto fc Saio 
Jl999h . 



• The red solid and short-dashed lines are the distributions adopted in K98, which are 
calculated with Eq.[2] and the mass ranges at Z = O.OO^I- Because of the metallicity 
dependent stellar lifetimes, the resultant lifetime distribution i?ia,t slightly depends on 
metallicity (solid lines, Z = 0.002; short-dashed lines, Z = 0.02). The two compo- 
nents appear at (md,tia) = (~ 2Mq, 0.5 — 1.5 Gyr) and (~ IMq, 2 — 20 Gyr), which 
correspond to the two types of binary systems, the MS+WD and RG+WD systems, 
respectively. A kind of bi-modality, which is observationally proposed, has been natu- 
rally realized in our model. Note that the spiked around ~ 1 Gyr are due to the change 

^The same as in Fig.l of KOO, but multiplied by the number fraction of WDs, / -^(f>{mi)dmi = 0.019. 
^In Fig.l of KOO, those are not seen due to a coarse mass grid. 
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in the stellar lifetime across the upper mass limit (~ 2Mq depending of metallicity) 
for the occurrence of the He core flash. 



The g reen dot-dashed lines are for the DD scenario, where the iTutukov fc Yungelson 



(Il994j )'s function is adopted in our chemical evolution model. The corresponding mass 
distribution (converted with the stellar lifetime function for Z = 0.02) spans a range 
of 1 — SmJ^I. Because is larger for larger mj, a typical lifetime is tia ~ 0.1 Gyr, 
the fraction with > 1 Gyr is very small, and the fraction with > 10 Gyr is almost zero. 

The blue long-dashed lines are for the MROl-like model, calculated with Eq. P, but 
with the same IMF as in our model. The mass distribution (converted with the stellar 
lifetime function for Z = 0.02) extends over — 8Mq having a peak at ~ 1.5Mq. 
The decrease of the rate toward larger m2 is due to the IMF, and the decrease toward 
smaller m2 is caused by mi + m2 > 3Mq. A typical lifetime is tja ~ 0.3 Gyr, which 
is longer than the DD model. The fraction with > 1 Gyr is two times larger than the 
DD model, and the fraction with > 10 Gyr is small but not zero. 



For each SN la progenitor model, the distribution function of the companion's mass and 
lifetime, i.e., the SN la rate in IMq of simple stellar population is given in Tables [T]l3l In 
our SD scenario, the SN la rate is set to be at [Fe/H] < —1.1. 



•^The drop of i?ia,m in the DD model is caused by the same reason as the spikes. 

Table 4. Upper and lower mass of the binary companion of SN la progenitor systems 
depending on the iron abundance, in the mass unit of Mq. 



[Fe/H] 


-1.1 


-1.0 


-0.7 





0.4 




0.9 


0.9 


0.9 


0.9 


0.8 




0.9 


1.5 


2.0 


3.0 


3.5 




1.8 


1.8 


1.8 


1.8 


1.8 




1.8 


2.6 


4.0 


5.5 


6.0 
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2.2. A New SN la Model 

Hachisu et al. (2008, hereafter HKN08) introduced the stripping effect, where a part of 
the envelope mass of the companion stars is stripped by the interaction with the WD winds. 
The efficiency is described by one parameter c, which depends mainly on the wind velocity. 
The effect of geometry can be given as a function of the mass ratio, and is included in c. 
The value of c can be determined from the observed binary systems such as supersoft X-ray 



sources, and is ~ 1.5 — 10 for one source (IHachisu fc Katdl2003abl ) and will be ~ 3 on the 
average. 

Although the binary evolution is calculated only for Z = 0.02 in HKN08, the wind 
velocity depends on the metallicity, and the case with higher metallicity should be similar 
to the case with larger c (Hachisu 2007, private communication). For Z = 0.002, the mass 
stripping is supposed to be inefficient, and thus the mass ranges of companion stars are the 
same as in K98 and KOO. For higher metallicity, the mass ranges are wider. Here we take 
the results with c = 0, 1, and 3 for Z = 0.002, 0.004, and 0.02, respectively. In the chemical 
evolution model, the upper and lower mass of the binary systems, and m^, are given as 
a function of [Fe/H] (Table H]), since iron is the most effective element for the opacity of the 
WD winds. While the metallicity effect is introduced as the simple cut-off the SN la rate in 
K98 and KOO, the metallicity effect is included in the parameter space of the binary systems 
that leads to SNe la in the present work. 

Figure [2] shows the SN la rate in the simple stellar population as functions of compan- 
ion's mass TTid (upper panel) and lifetime tja (lower panel) for the metallicity Z = 0.002 
(blue), 0.004 (green), 0.02 (orange), and 0.05 (red). The sohd and dashed lines indicate the 
MS-I-WD and RG+WD systems that correspond to the young and old population of SNe la, 
respectively. The total number of SNe la is determined to reproduce the chemical evolution 
in the solar neighborhood ([6rg, &ms] = [0.023, 0.023] at Z = 0.004, see §3 for the details), 
and also adopted for the other types of galaxies (§4). We should note that this function 
is different from HKN08 because the adopted IMF, stellar lifetimes, and binary parameters 
are different. This function in Figure [2] is better for chemical evolution models since it is 
calibrated to meet the observations in the solar neighborhood. 

Figures [3] show the summation of the young and old components of the SN la rate 
as functions of companion's mass (upper panel) and lifetime (lower panel). The bimodal 
distribution is realized at (md,tia) = (~ 2Mq, 0.1 — 0.5 Gyr) and (~ IM©, 1 — 3 Gyr) 
corresponding to the MS+WD and RG+WD systems, respectively. The lifetimes at the two 
peaks depend on the metalhcity: tja ~ 0.5 and 3 Gyr for Z = 0.002 (blue short-dashed line), 
tia ~ 0.2 and 1 Gyr for Z = 0.004 (green solid lines), tja ~ 0.1 and 1 Gyr for Z = 0.02 (orange 
dotted line) and 0.05 (red long-dashed line). Note that the spikes appears around ~ 1 Gyr 
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because the stellar lifetime is substantially longer for smaller mass stars that undergo the 
He core flash. This enhanced rate might cause the SN la driven galactic winds in elliptical 
galaxies. 

At higher metallicity, the lifetime distribution extends to both shorter and longer life- 
times. Such a short lifetime as tia ~ 0.1 Gyr is caused by the stripping effect that enables 
the accreting WDs with massive (~ 5 — QMq) companions to explode. The lifetime as long 
as tia ~ 20 Gyr is caused by the elongation of the lifetime of small-mass companion stars, 
which is dominant in present-day elliptical galaxies. 

HKN08 model predicts the existence of various types of circumstellar matter (GSM) 
around the WDs just before making SN la explosions depending on the mass accretion rates 
(and thus on the binary parameters) at the last stages. (1) case WIND: If the accretion 
rate is high enough for the optically thick WD wind continue to blow, the stripped-off 
material forms GSM very near the WD because the velocity of GSM is rather low; this case 
corresponds to SNe 2002ic (e.g., Deng et al. 2003) and 2005gj (Aldering et al. 2006), where 
SN la undergoes circumstellar interaction. (2) If the accretion rate has decreased below the 
critical rate for the wind to occur at the explosion, the stripped-off material forms GSM with 
various speed, but it has been dispersed far from the WDs. Depending on the distance from 
the WDs and geometry of the materials, GSM features may or may not be observed as in SN 
2006X (Patat et al. 2007) or SN 2007af (Simon et al. 2007). We also note that the speed of 
the stripped-off materials is as low as 10 — 100 km s~^, so that the sweeping effect of very 
fast wind discussed by Badenes et al. (2007) is not expected. 



2.3. Metallicity Effect 



After our theoretical prediction of the metallicity effect, the SNe la in low-luminous 
galaxies have been paid attention. KOO commented on SN1895B and SN1972E in NGG 
5253, and SN1937G in IG 4182 that the metallicities of the host galaxies are higher than 



[Fe/H] = -1.1 
< -1 (Fig. 1 in 
galaxy with Mr 



Hamuy et al. 



Hamuv et al. 



2000 ) showed that no host galaxies are found at [Fe/H] 
200lh N1999aw is in the 

. 12.2 (jStrolger et al.ll2002l ). which roughly corresponds to [0/H] ~ —1.2 

jTremonti et allbooi) SN2005cg is in the host galaxy with M, = -16.75 and [0/H] = -0.4 



SN2007bk is in the galax y with [0/H] 



-0.59, but located at 

a distance of 8.7 kpc (iPrieto. Stanek. fc Beacornll2008l ). where [0/H] is estimated as ~ —0.7 



(see below for the details). iBadenes et al.l (120091 ) showed the stellar population near the su- 
pernova remnants (SNRs) of SNe la in the Small Magellanic Gould from the color-magnitude 
diagram. The mean stellar metallicities are [Fe/H] = —0.96, —0.94, —1.15, and —0.80, but 
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there is also a young metal-rich population as they commented. From the statistical studies 
of a large sample such as the Sloan Digital Sky Survey and th e Supernova Legacy Survey, 



Gallagher et al. 



2005, 2008 



no such met alii city effect has been seen in observations yet (e.g. 

Howell et al.lboool ). and rather op posite effect (i.e., decreasing of S N la rates toward higher 



metallicity) has been suggested by ICooper. Newman fc YanI (120091 ). 



Note, however, that the following uncertainties are involved in these arguments; i) The 
observation is the emission weighted gas phase [0/H], and can be different from the [Fe/H] of 
the progenitor system; gas metallicity can be larger then stellar metallicity in the chemically 
evolved and closed system, or can be smaller due to the dilution by gas accretion and/or inho- 
mogeneous effects (i.e., stars tend to be formed from metal-enriched g as), ii) The metallicit y 
correlates with the galaxy mass (the mass-metallicity relation, e.g., iTremonti et al.l 120041 ). 



but the scatter is larger for smaller galaxies, iii) The mean stellar [a/Fe] correlates with the 
gala xy mass in the sense that ] a/Fe] is larger for more massive galaxies ([a/Fe] ~ — 0.3, 
e.g.. lThomas fc Marastonll2003l ). although the errorbar is quite large, iv) There exits a radial 
gradi ent of metallicity, where the metallicity is lower in outside (e.g., iKobayashi fc Arimoto 
19991 ). v) There exits a metallicity distribution of stars even at a certain location in galaxies 

(Fig. SD. 



Foerster fc Schawinski I (120081 ) reported that there is no radial dependence of SN la rate 



found in early-type galaxies. However, this work also did not reach the metallicity threshold. 
The metallicity gradient of early- type galaxies is typically A[ Fe/H]/Alogr ~ —0.3, and th e 
metallicity at the effective radius is < [Fe/H] >c~ —0.3 to (IKobayashi fc Arimotdll999l ). 
The metallicity effect would be seen at r > lOre in intermediate-mass galaxies. In dwarf 
galaxies, SNe la can be found at r ~ 50 kpc, because the effective radius does not become 



so small ( Binggeli et al. 



galaxies (iSpolaor et al 

[0/H] ~ [Fe/H] ~ -0.7 with A [Fe/H] /A log r/n 



198J), and because the metallicity gradient is shallower for smaller 
2OO9I ). For SN2007bk, the local metallicity is roughly estimated as 

1.5 - 2 kpc. 



-0.2 and ~ 



3. [X/Fe]-[Fe/H] Relations 

Using the one-zone chemical evolution model, we show the evolution of elemental abun- 
dance ratios in the solar neighborhood for different SN la models. The star formation history 
is determined from the following observations in the solar neighborhood. The formulation 
is described in KOO; we use the model that allows the infall of primordial gas from outside 
the disk region. The SFR is assumed to be proportional to the gas fraction. For the infall 
rate, we adopt the same formula as in K98 and K06 that is proportional to texp[— ;^]. The 
input parameters are the Galactic age of 13 Gyr and infall timescale Xi = 5 Gyr. The star 
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formation timescale Ts = 2.2 Gyr is determined from the present gas fraction /g = 0.15 (see 
KOO for the definition). The same parameters are set for the other SN la models, because 
the contribution of SNe la on the gas fraction is neghgible. The nucleosynthesis yields as 
functions of mass, metallicity, and explosion energy are given in K06 for HNe, SNe II, and 
SNe la. Photometric evolution is also calculated as in KOO. 

Figure H] shows (a) the SFR, (b) the age-met allicity relations, and (c) the metallicity 
distribution function (MDF) for different SN la models. The resultant SFR peaks at t ~ 8 
Gyr, but does not very much change for t ~ 5 — 13 Gyrs, and the present SFR is consistent 
with the observational estimate. The red solid, cyan dotted, green short-dashed, blue long- 
dashed, and magenta dot-dashed lines are for our new model with and without the metallicity 
effect, the DD, MROl-like, and K98 models, respectively. All models can give excellent 
agreement with these observations. Compared with the K98 model (magenta dot-dashed 
line), our new SN la model (red solid line) gives slightly slower increase in [Fe/H] after 
[Fe/H] ~ — 1 (panel b), and slightly larger number of metal-poor stars around [Fe/H] ~ — 1 
(panel c). This is because the metallicity effect i s included not with a sh arp cut-off but with 



a gradual decrease in our new model . Recently. iNordstrom et al.l (20041) showed a n arrower 



MDF than the plotted observations (lEdvardsson et al.lll993l : IWyse fc Gilmord 119951 ) . which 



may suggest that a sharp cut-off exist in the SN la rate at [Fe/H] ~ — 1. 

The total number of SNe la is chosen to meet the MDF; For the DD model, the total 
fraction of SNe la is set to be 0.001. For the MROl-like model, the binary fraction (see §5 for 
the definition) of A = 0.035 is adopted. In our models, the binary parameter b denotes the 
fraction of primary stars that eventually explode as SNe la for the RG+WD and MS+WD 
systems. The combination of 6rg and 6ms can be determined from the metal-rich end of 
the MDF and the [0/Fe]-[Fe/H] relation at [Fe/H] > - 1. For our new SN la model, b is 
normalized at Z = 0.004, and [&rg, &ms] = [0-023, 0.023] are chosen for the best model. (Note 
that the models with [0.015, 0.028], [0.010, 0.032] and [0.005, 0.035] are also possible.) Since 
the mass range expands for higher metallicity, the fraction of SNe la is 0.040 and 0.032 at 
Z = 0.05 for the RG+WD and MS+WD systems, respectively. For the case without the 
metalhcity effect, the companion's mass ranges for Z = 0.02 and [6rg, &ms] = [0-027, 0.027] 
are adopted. For the K98 model, [b^c, 6ms] = [0-02, 0-04] are chosen in K06 for the updated 
nucleosynthesis yields- 

In the solar neighborhood, elemental abundance ratios are observed for a number of 
stars over a wide range of metallicity, which gives one of the most stringent constraints on 
the progenitors of SNe la- Figures 013 show the evolutions of elemental abundance ratios 
[X/Fe] for O, Mg, Mn, and Zn, against the iron abundance [Fe/H]. 

In the early stage of the galaxy formation, only SNe II explode, and [a/Fe] stays con- 
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stant. Because of the delayed Fe production by SNe la, [a/Fe] decreases toward 0. The 
observation s with UV OH hues show a monotonic in crease in [0/H] towards the lower 
metaUicity (llsraelian et al.l Il998l : iBoesgaard et al.lll999l ). but these are largely affected by 
the 3D effect of the atmosphere model and are not plotted here. The O observations with 
the forbidden line and infrared lines are more reliable, and show the plateau rather than the 
increase as well as the other a-elements (Mg, Si, Ca, and Ti). 

Manganese shows a characteri stic behavior. Observational data show that [Mn/Fe] 
increases toward higher metaUicity. iFeltzing. Fohlman. &: Bensbyl ( 120071 ) concluded that the 
Mn trend is likely to be due to the metaUicity dependent SNe II, but the trend can be 
naturally reproduced by the delayed enrichment by SNe la since SNe la produce [Mn/Fe] 
> (K06). This is consistent with that this increasing trend of [Mn/Fe] appears from the 
same [Fe/H] as the decreasing trends of [a/Fe]. In the observation, a larger increase of 
[Mn/Fe] may be seen at [Fe/H] ~ than our models, which may be due to the metaUicity 
dependence of SN la yields, which is not included in our models. 

The species of zinc depend on the metaUicity. At low metaUicity, ^^Zn is produced 
by complete Si- burning in HNe (jUmeda fc Nomotdl2002l ). which are assumed to be a half 
of massive (> 2OM0) SNe II with the lifetime of < lO'^ yr (K06). At high metaUicity, 
neutron-rich isotopes of zinc ^^~™Zn are produced by neutron-capture in He and C burning, 
which is larger for higher metaUicity (Fig. 5 in K06) and ejected by SNe II. Because of the 
combination of the lifetime and metaUicity effects, [Zn/Fe] shows an interesting track. 



• In the DD model (green short-dashed line), [a/Fe] decreases and [Mn/Fe] increases 
too early and too quickly from [Fe/H] ~ —2 compared with the observations. This is 
because the lifetime of the majority of SNe la is shorter than 0.3 Gyr (Fig. [T]). The SN 
la lifetime is so short that [Zn/Fe] decreases by 0.25 dex from [Fe/H] ~ — 2 to ~ —0.6. 
Because the SN la rate sharply decreases for longer Ufetime, and because the SFR does 
not change very much, the contribution of HNe becomes larger than SNe la at later 
time. This results in the [Zn/Fe] increase from [Fe/H] ~ —0.6 to ~ 0, which is not 
seen in the observations. 

• In the MROl-like model reproduced with our code (blue long-dashed line), the decrease 
in [a/Fe] and the increase in [Mn/Fe] are seen at [Fe/H] ~ —2 as early as in the DD 
model. The slope of [a/Fe] and [Mn/Fe] against [Fe/H] is a little shallower than the 
DD model because the fraction of lifetime with > 1 Gyr is larger. The [Zn/Fe] decreases 
only by 0.15 dex from [Fe/H] ~ —2 to ~ —0.6. 



Our SN la model with the metaUicity effect (red solid line) gives the best agreement 
with the observations. The decrease in [a/Fe] and the increase in [Mn/Fe] are deter- 
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mined not from the lifetime effect but from the metaUicity effect on the SN la rate. 
Since the lowest metallicity to produce SNe la is [Fe/H] = —1.1, the evolutionary 
change is caused at [Fe/H] ~ —1 by the companion stars with Mms,u = 2.6Mq after 
the main-sequence lifetime of ~ 0.5 Gyr. In the K98 model, the metallicity effect is 
included more sharply, and thus the evolutionary change appears a bit more sharply. 

Our model predicts an interesting evolution of [Zn/Fe] around [Fe/H] ~ —1, as a 
contrast to other models; [Zn/Fe] increases from [Fe/H] ~ —1.5 by the metallicity 
effect on Zn yield, and decreases from [Fe/H] ~ — 1 by the contribution of SNe la. 
This t rend seems to be co nsistent with the obse rvations, although the scatter is not 



small fiNissen et al.l (120041 ) at [Fe/H] < -1 and iNissen et al.l (120071 ) with the NLTE 
correction at [Fe/H] > —1). Further observations of these elements at [Fe/H] > — 1.5 
is critically important in order to identify the SN la progenitors and to clarify the 
metallicity dependence on the nucleosynthesis yields. 

If we do not include the metallicity effect (dotted line in Fig. E]), [a/Fe] starts to 
decrease at [Fe/H] ~ —2. This corresponds to the shortest SN la lifetime of ~ 0.1 Gyr 
for Mms,u = 5.5M0. This is too early to be compatible with the observations. In the 
K98 model, the [a/Fe] decrease starts at [Fe/H] ~ — 1.5 (see Fig. 3 in K98), because 
the shortest lifetime (~ 0.5 Gyr for Mms,u = 2.6Mq) is longer than our new model. 
In other words, without the metallicity effect, the existence of the young population 
of SNe la cause the [a/Fe] decrease at much smaller [Fe/H], and the discrepancy with 
the observations becomes larger. 

If there is the metallicity effect, these results do not depend very much on the c pa- 
rameter (Fig. [7]). Even with c = 10 that leads to the shortest lifetime of 0.03 Gyr, 
the difference is not visible. This is because the majority of stars have low metallicity 
a.s Z < Zq, and the metal- rich component can contribute only after ~ 10 Gyr. If the 
stripping is effective even for Z = 0.002 (long dashed line), the results does not change 
very much. Therefore, the evolutionary change in the [X/Fe]-[Fe/H] relations is simply 
caused by the limit of the WD winds, as introduced in K98. 



Figure [8] shows the SN la rate history in the solar neighborhood. Even with the same 
present SN la rate, different SN la models produce different SN la rate histories. For the 
DD model, the lifetimes of majority of SNe la lifetime are short, and the SN la rate shows a 
rapid increase at t < 1 Gyr and then a constant until the present time. For the SD model, 
the lifetimes are a bit longer, and the SN la rate shows a slower increase at t < 5 Gyr. 
With the metallicity effect, SNe la start to occur at 3 Gyr when [Fe/H] reaches -1.1 (FigJDD). 
Compared with our previous work, our new model shows a more gradual increase due to the 
smoother SN la lifetime distribution as a function of metallicity. 
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In Figure [9], the supernova rate is plotted against the metalhcity indicators: (a) the iron 
abundance, (b) the gas-phase metalhcity, and (c) the mean stellar metallicity (see Eq.[15] in 
KOO for the definition). Since the SN la rate is assumed to be at [Fe/H] < —1.1, the rapid 
decrease is seen at [Fe/H] ~ —1 after the short time-delay of the chemical enrichment (panel 
a). Because of the oxygen enhancement at lower metallicity, the decrease is seen at higher 
[0/H] as ~ —0.6. In dwarf galaxies, however, the decrease will be seen at [0/H] ~ —1 since 
the observed dwarf galaxies have low [0/Fe]. For the mean metallicity (panel c), which is 
observed with the absorption lines of integrated stellar populations, the dec rease is seen at 



lower < [Fe/H] > as ~ —1.2. With the mass- metallicity relation of galaxies (ITremonti et al. 



20041 ). the metallicity threshold corresponds to Mb = ~15 mag and 3 x lO^M, 



0- 



Mannucci et al.l (120061 ) suggested that ~ 50% of SNe la sho uld be composed of the 



systems with lifetimes as short as 0.1 Gyr. On the other hand, iMatteucci et al. J2006h 
claimed that the fraction should be less than ~ 35% in order to explain the [0/Fe]-[Fe/H] 
relation and iron abundance in the intracluster medium. In our model, although 50% of SNe 
la eventually come from the MS+WD systems, the lifetime shorter than 0.1 Gyr is possible 
only for such higher metallicity as Z > Zq. Figure [TO] shows the time evolution of the 
cumulative function of SNe la lifetime in the solar neighborhood, where the star formation 
and chemical enrichment histories are taken into account. At the galactic age of t = 3 Gyr, 
only the MS+WD system can produce SNe la, and thus the lifetime of all SNe la is less 
than 1 Gyr. The RG+WD systems contribute at later time. At present (t = 13 Gyr), the 
fraction of SNe la that have shorter lifetime than 1 Gyr is ~ 50%. The total fraction with 
< 0.1 Gyr is only 10% in our model. 



Table 5. The input parameters of galaxy models: timescales of the star formation and 

infiow in Gyr. 



r, [Gyr] [Gyr] 



E 


0.1 


0.1 


SOa/Sa 


1.8 


2.0 


Sab/Sb 


1.75 


4.5 


Sbc/Sc 


2.0 


8.3 


Scd/Sd 


2.1 


23.0 
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4. Galactic Supernovae Rate 

4.1. SN la rate in Ellipticals 

The present SN la rate in elliptical galaxies gives also a stringent constraint on the SN 
la progenitor models. Figure [11] shows the SN la rate history in ellipticals, which is similar 
to Fig. 5 of KOO, but for the cosmological parameters of Hq = 70 km s~^ Mpc~^, Qq = 0.3, 
Ao = 0.7, the galactic wind epoch of 1.5 Gyr, and the galactic age of 13 Gyr. (Note that 
the observed SN la rate is proportional to Hq.) The star formation history is assumed as 
in KOO; the bulk of stars in ellipticals are formed at z >3, and thus having ages older than 
10 Gyr. The adopted inf all and star formation tim escales are summarized in Table [51 The 



present B-V color of 0.9 ([Roberts fc Hayneslll994l ). the mean stell ar metallicity of ~ ~0-2 
( Kobayashi fc Arimoto 19991 ). and the mean stellar [a/Fe] of 0.2 ( Spolaor et al. 2009) are 
reproduced. The present stellar mass-to-light ratio is M/Lk = 1.8 and M/Lb = 9.3 in our 
model. The resultant SFR is similar to the SN 11 rate (black dot-dashed line). 

• In the DD model (green short-dashed line), th e SN la lifetimes a.re too short to explain 
the observed SN la rate at the present epoch f Cappellaro et al.lll999l ). Toward z ~ 4, 



the SN la rate rapidly increases, and the SN la rate at 2; ~ 4 is 10 times larger than 
the present rate. 

In the MROl-like model (blue long-dashed line), the fraction of SNe la with tja > 10 
Gyr is larger than the DD model, so that the present rate is not too small compared 
with the observation. Toward 2 ~ 3, the SN la rate gradually increases, and the SN 
la rate at 2; ~ 1 and ~ 3 is larger by a factor of 2 and 4, respectively, than the present 
rate. 

In ellipticals, the chemical enrichment takes place so early that the metallicity effect 
on SNe la is not effective. Thus the SN la rate depends almost only on the lifetime. 
Our model (red solid line) includes the RG-fWD systems with tja > 10 Gyr as the SN 
la progenitors, and thus the observed SN la rate in ellipticals can be well reproduced. 
Toward higher redshift, the predicted SN la rate does not change very much (in the 
supernova units, SNu). 

In the K98 model (see Fig.5 in KOO), the lifetime distribution shows a gap between the 
two components of the MS-I-WD and RG+WD systems (Fig. [1]), so that two peaks 
are clearly seen in the SN la rate in ellipticals. Such peaks are not clearly seen in our 
new SN la model because the lifetime ranges get gradually wider for higher metallicity 

(Fig. m- 
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The sharp peak at t ~ 1.5 Gyr is caused both by the initial star burst assumed in the 
SFR and by the spikes seen in Figure [3l which are caused by the enhancement of the stellar 
lifetime from the occurrence of the He core flash. The galactic winds might be driven by 
this enhancement of the SN la rate in elliptical galaxies. Therefore, we set the galactic wind 
epoch of 1.5 Gyr in our model of elliptical galaxies. 



4.2. SN la rate depending on Galaxy Type 

Since the star formation history is different for different types of galaxies, the SN la 
and II rates depend on the morphological type of galaxies. In the observations, the present 
rates are available for various types of galaxies, which can give a constraint on the SN 
la models. In Fig. 4 of KOO, the observed ratios between SN la and II rates have been 
reproduced with our SN la raodel. We show the similar figure to compare with the updated 



observation (IMannucci et al.l 120051 ). We should note, however, that the absolute rates per 
B-band luminosity involves an uncertainty from dust extinction, and the rates per mass 
involve an uncertainty in the mass-to-light ratio. 

We update the evolution models of four types of spiral galaxies, in order to meet the 
observational constraints on the present colors and the gas fractions, assuming the galactic 
age of 13 Gyr. The adopted infall and star formation timescales are summarized in Table 
[51 Figure [12] shows the time evolution of (a) the SFR, (c) the gas fraction per luminosity, 
and (d) the B-V color for four types of spirals (red dotted line for SOa/Sa; yellow solid 
line for Sab/Sb; green long-dashed line for Sbc/Sc; blue short-dashed line for Scd/Sd). In 
earlier-type of spirals, star formation takes place earlier, and thus the earlier-type of spirals 
have older ages of stellar populations, smaller gas fractions, and redder colors at present. In 
the panel (b), the upper and lower four lines respectively show th e SN II and la ra t es per 



B-band luminosity (in SNu). The observational data is taken from ICappellaro et al.l (119991 ) 
for SOa-Sb and Sbc-Sd. The present SN Il-I-Ibc rate is larger for later-type of spirals, while 
the SN la rate in SNu is not different so much among the various types of galaxies, which 
are roughly consistent with the observation. 

Figure [13] shows the present supernova rates per K-band luminosity (upper panel) and 
per mass (lower panel) against the morphological type of galaxies. The supernova rates are 
larger for later-type of spirals because the SFR is higher. Almost parallel relations between 
the observed rate and the galaxy type are seen for SNe II and Ibc, which may suggest 
the binary fraction is universal being independent of the galaxy type. Among core-collapse 
supernovae, we assume that 85% and 15% are SNe II (blue dashed line) and SNe Ibc (green 
dotted line), respectively. Thus, it seems to be reasonable to adopt the same b parameters as 
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SNe la for all types of galaxies. As a result, our SN la model (red solid line) gives excellent 
agreement with the observed SN la rates. In our SN la models, there are two types of 
progenitor systems; The SN la rate is larger for later-type spirals, which is due to the young 
population of the MS+WD systems. The slope of the SN la rate against the galaxy type 
is flatter than that of the SN II and Ibc rates, which is due to the old population of the 
RG+WD systems. 



The observational data is taken from iMannucci et al.l (120051 ). where the adopted mass 



to-light ratios are M/Lk ~ 0.8,0.7,0.5, and 0.4 for E/SO, SOa/b, Sbc/d, and Irr in their 
Table 2. In our models, the "stellar" mass-to-light ratios are M/Lk ~ 1.29,1.14,1.07, and 
1.03 for SOa/Sa, Sab/Sb, Sbc/Sc, and Scd/Sd, which are ~ 2 times larger than their ratios. 
This difference might come from the difference in the IMF. We thus plot the observational 
data multiplied by a factor of 2 for the rates per mass. 

Figure [T3] compares the galactic supernova rates for different SN la models. 



• In the DD model (green short-dashed line), the lifetime of the majority of SNe la is 
short, no large difference between the SN II and la rates is seen in the trend against 
the galaxy type, which is not consistent with the observations. 

• The MROl-like model (blue long-dashed line) and our SD models (red solid line and 
cyan dotted line) give almost the same trends, regardless of the metallicity effect, as 
far as the lifetime distributions in this work and K98 are applied. 



4.3. Cosmic Supernovae Rate 



Galaxies have various timescales for star formation and chemical enrichment, and the 
occurrence of SNe la depends on the metallicity therein. Therefore, we should calculate the 
cosmic supernova rate by summing up the supernova rates in spirals and ellipticals with 
the ratios of the relative mass contribution in the Universe (see KOO for the detail). As in 
Fig. 6 of KOO, the resultant cosmic SFRs shows a peak at ^ ~ 2. Figure [TBI shows the cosmic 
supernova rates with our new SN la model (solid lines), as a composite of those in spirals 
(blue lines) and ellipticals (red lines) in SNu (upper panel) and in the unit of the comoving 
density (lower panel), respectively. The SN II rate (green dot-dashed lines) traces the star 
formation history, and shows a pea.k at z ~ 2 in the comoving d ensity, being consistent with 
the observations (ICappellaro et al.ll2005uBotticella et al.ll2008l ). 



As in spirals, the total SN la rate (green solid lines) gradually decreases toward z ~ 2.6 
in the unit per luminosity, while the total SN la rate slightly increases to 2; ~ 1 and then 
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decreases in the unit per volume. In ellipticals, the chemical enrichment takes place so 
quickly that the metallicity is large en ough to produce SNe la at z > 2. These are broadly 



consistei it with the obsery ations (e.g., iPain et al.ll2002l . Supernova Cos mology Project, 38 



SNe la; iTonry et al.l |2003| . High-z Supernova Search T eam, 8 SNe la; iDahlen et al.l 12004 



Great Observatories Origins Deep Sur vey, 25 SNe la: iBarris fc Tonryl |2006| . Institute for 
Astronomy Deep Survey, 98 SNe la?; iNeil et al.l l2006l . Supernova Legacy Survey, 58 SNe 
la). Around z ~ 0.5, t he total SN la rate may be smaller than the observations, although 
Poznanski et al.l (120071 ) shows a significant lower rate with the Subaru Deep Field. 



The cosmic SN la rate history for different SN la models are also shown in Figures [T5] 



and 



• If we do not include the metallicity effect (dotted lines in Fig. [15]), the SN la rate in 
spirals does not decrease so much, and the total SN la rate per luminosity is almost 
constant. The rate per volume slightly increases toward higher redshifts, and is larger 
than the present rate by a factor of ~ 2 at z ~ 2. In the observations, the decrease in 
the total rate is seen, which supports the existence of the metallicity effect. 

• In the DD (short-dashed lines in Fig. [T6l) and MROl-like (long-dashed lines in Fig. 
[T6|) models, the lifetime of the majority of SNe la is so short that the rate in ellipticals 
keeps on increasing at z < 4, and that the rate in spirals shows a peak at 2 ~ 1.5. As a 
result, the total SN la rate per luminosity is almost constant, and the rate per volume 
keeps on increasing toward higher redshifts. The total SN la rate is larger than the 
present rate by a factor of ~ 2 and ~ 4 at 2; ~ 1 and 2, respectively. 

• In the K98 model (see Fig. 6 in KOO), two peaks are seen in the SN la rate history 
in ellipticals, which are not clearly seen in our new SN la model, because the lifetime 
ranges get gradually wider for higher metallicity (Fig. [3]). 



The observed cosmic SFR shows an increase by a factor of ~ 10 from 2; = to ~ 1 



[e.g., iSchiminovich et al.ll2005l ). and such a large increase is not seen in the composite SFRs 
of these models for ellipticals and spirals. Such a large increase in the SFR results in the 
rapid increase in SN II and la rates from 2; = to ~ 1 that may be seen in the observed SN 
la rates. By modifying the SFRs in spirals as well as including the contribution of irregular 
galaxies, there can be a better fit (ICalura &: Matteuccil 120061 ). However, we don't apply 
such parameter-tuning in this paper because other effects are more important in the cosmic 
chemical enrichment: i) the number evolution of each type of galaxies, ii) merging of galaxies, 
iii) internal structure, namely, metallicity gradients within galaxies, iv) the contribution of 
dwarf galaxies. We will rather show the cosmic supernova rates with fully self-consistent 
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cosmological simulations (iKobayashi et al.l 120071 : lKobayashill2008l ). Nonetheless, our result 
with the simple assumptions is robust; the SN la rate is high in ellipticals and drops in 
spirals at high-redshift. 



5. Discussion 

5.1. SN la Formulation 

In galactic chemical evolution models, the following two formulations have been pro- 
posed. Their parameterizations are very different, and it would be worth summering here. 
Observations of binary systems could give a clue to improve these formulations. 

(1) The formulation of the SN la rate has been first proposed by Greggio & Renzini 
(1983). The SN la rate is described as 

7^I, = A / (piMs) f{fx)^it-rM,)diidMB. (1) 

The first integral is for the total mass of the binary Mb = Mi + M2, and the second integral 
is for the mass fraction of the secondary star, /i = M2/MB. In MROl, the Scalo IMF 
is adopted, and the distribution function of the mass fraction is described as /(/i) = 
2^+^(1 + 7)/i'*'. 7 = 2 is adopted in their successful models, which gives weight for massive 
secondaries in the distribution function of mass ratios (see Fig. [T^ . The upper and lower 
limits of integrals are MB,inf = max[2M2(t), Mem] and Mb,sup = ^bm/2 + M2(t), where 
the limits to the total mass of the binary are M-Bm = SM© and Mbm = IQMq. Therefore, 
the initial mass range of the primary stars is 1.5Mq < Mi < 8Mq and the secondary 
stars are AGB stars with the initial mass of M2 < 8Mq. The lower limit of the primary, 
however, should not l.SM© but ~ SMq, since the primary should be a C+0 WD. Although 
A is called as the binary fraction, A a ctually indicates the frac tion of binary systems with 



suitable separation to become SNe la (iGreggio fc Renzinilll983l ). 



(2) Based on the SD scenario, an alternative formulation has been proposed by K98 and 
KOO. The SN la rate is given as 



I — 0(m) dm / — 'ijj{t — Tm) (pdi'm) dm. (2) 

max [jTtp ^, mt] ^ J ma,x[m^ i,mt] ^ 



These integrals are calculated separately for the primary and secondary stars (mp and ma) 
respectively. In our works, the Salpeter IMF is adopted, and the distribution function of 
the companion's mass is assumed to be the power-law as 0d(^) m~°-^^. The slope —0.35 is 
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chosen to meet the observations (e.g.. iDuquennoy &: Mayorlll99ll ). and such negative slope 



gives the distribution function weighted for less-massive secondaries (see Fig. [19]). 0d is 
normalized to unity in the integrated mass ranges as J^^'^ ^4>d{'^)dm = 1. Since the 
primary stars are C+0 WDs, the initial mass range of the primary stars is rrip = 3 — 8Mq. 
The initial mass ranges of the companion stars, rrid/ and md,^, are given by the simulation 
of binary evolution (Hachisu et al. 1999ab) for the RG+WD and MS+WD systems. The 
binary parameter b denotes the total fraction of primary stars that eventually explode as 
SNe la, and is determined from the chemical evolution in the solar neighborhood, b includes 
the binary fraction and the efficiency for each binary system. In addition, the metallicity 
effect is taken into account (§2). The lowest metallicity to produce SNe la was set to be 
[Fe/H]= —1.1 both in K98 and KOC0. In our new model (§2.2), the mass stripping effect is 
taken into account, and and md,u are given as a function of metallicity as in Table HI 

Strictly speaking, nid/ and md,u depend on mi, but it can be neglected. For Mwd = 
0.7Mq, the mass range is a little narrower, and thus the shortest lifetime is 0.2 Gyr that 
is longer than 0.1 Gyr for Mwd = l.OM© {Z = 0.02). The range of the orbital period is 
much narrower, so that the contribution of Mwd < l.OM© is smaller. On the other hand, 
for Mwd = I.IM©, the shortest lifetime is a little shorter, but the number of such WD 
supposed to be very small. As a result, the [X/Fe]-[Fe/H] relation is mainly determined by 
the contribution for Mwd = l.OM©. 

The results slightly depend to some extent on the stellar lifetime function of low and 
intermediate mass stars. Figure [T7] shows the adopted stellar lifetime function as a function 
of initial mass and metallicity, which is provided by Kodama & Arimoto (1997). The stellar 
evolution tracks are calculated with the code described in Iwamoto & Saio (1999) with the 
core overshooting. Because of the coarse mass grid, there is a tiny irregularity at l.SM© 
(Z = 0.004,0.008), 2.OM0 {Z = 0.02), 2.2M0 (Z = 0.0001-0.002), and S.SM© {Z = 0.0002) 
in the table, which we have smoothed wit h interpolation. The stellar li fetime function overall 



agrees with more recent models such as iKarakas fc Lattanzid (120071 ). The long and short 



dashed-lines are for the analytical functions adopted in MROl and HKN08, respectively, 
where the metallicity effect is neglected. With those analytical functions, the stellar lifetime 
is systematically shorter. 



We adopt the Salpeter IMF with a slope of —1.35. iKroupal (120081 ) showed a similar 
slope (—1.3 ± 0.5) for m > O.5M0, which is flatter than in the Scalo IMF. For m < O.5M0, 
the slope is shallower and the number of brown dwarfs are smaller than in the Salpeter IMF. 



''The note in MROl is incorrect. See K98, £.16 in the right column, p. 156, and KOO, £.15 in the right 
column, p. 27. 
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In chemical evolution models, however, metal enrichment is proceeded only by stars with 
m > 0.5Mq. If a suitable value is chosen for the lower mass-cut of the IMF, the Salpeter 
IMF can give equal results with the Kroupa IMF. Therefore, we adopt the single slope for 
a mass range of m = 0.07 — 50Mq. This does not necessarily mean that there is no star at 
m < O.OTMq. The requirement is that the mass fraction of m = 0.5 — SOM© is 45%. 



5.2. Binary Parameters 

In the two formulations, Eqs.[T] and [2], the treatment of binary systems is different 
for the following three points: i) binary fraction, ii) mass ratio, and iii) binary IMF. These 
differences may reflect the formation processes of binary systems. Theoretically, the prompt 
fragmentat ion scenario s eem to be more feasible than captures or delayed breakup of accre- 



tion disks (jTohlindl2002l ). and these scenarios will be tested with the observed distribution 
functions of the binary parameters. In chemical evolution models, however, these uncertain- 
ties can be cancelled out by the choice of b or A. 

i) Binary fraction — The total number of SNe la, namely, the binary fraction A or 
binary parameters b, could be determined from the binary population synthesis. However, 
the evolution of binary systems is complicated and the effect of the WD winds is not included. 
In Eq. [2], we treat b as parameters being determined for the RG+WD and MS+WD systems, 
independently, from the requirement to meet the observed MDF of G-dwarf stars and the 
[0/Fe]-[Fe/H] relation in the solar neighborhood. In K98, we adopted 6rg = ^ms = 0.04 from 
the Kolmogorov-Smirnov (KS) test of the [0/Fe]-[Fe/H] relation, but in KOO, we adopted 
[&RG,&Ms] = [0.02,0.05] from the testi. In K06, we adopted smaller values of [0.02, 0.04] 
because Fe are also produced by HNe in the updated nucleosynthesis yields. For our new SN 
la model, [6rg, &ms] = [0.023, 0.023] are chosen at Z = 0.004, which becomes [0.040, 0.032] 
at Z = 0.05, as discuss ed in S 3 . In t otal, the fraction of SNe la is 5 — 7%, which is consistent 



with the constraints in iMaoa (120081 ) except for the X-ray observations 



In Eq.pp, the definition of A is totally different from that of b, and A can be simply 
determined from the MDF. This is misunderstood in Matteucci & Recchi (2001), where 
they claimed that they failed to reproduce the K98 results (Mod.l in MROl). We would 
demonstrate in Figure [TSlthat we can reproduce their and our results. For the comparison, we 
adopt the same IMF, stellar lifetimes, and nucleosynthesis yields as in our models. Adopting 
MROl's formula Eq.[T] and parameters, we almost reproduce Mod.l and Mod. 3 in MROl 
(short- and long-dashed lines, respectively). Here we adopt A = 0.035 to meet the MDF, 



^There are two miss-prints in Appendix of KOO, not [6ms — 0.02, 6rg — 0.05] but [6rg — 0.02, bus = 0.05]. 



-22- 



which is better than A = 0.05 originally adopted in MROl. 



In Mod.3 in MROl, A = 0.02 and 0.05 are adopted for the RG+WD and MS+WD 
systems, respectively, which are the same values of the binary parameter b in KOO. However, 
since 0d is normalized to unity not for m = 0.07 — 5OM0 but for ma ~ 1 — 3Mq, the value 
of A should be ten times larger than the value of b. Adopting A = 0.45, we reproduce our 
original K98 model (solid line) even with the MROl formula (dotted line). With the MROl 
formula, the slope of [0/Fe] is slightly different, and [0/Fe] around [Fe/H] ~ —0.5 is slightly 
higher than the K98 model. The best choice of A depends on the two other assumptions 
(i and ii), and varies in the range of A = 0.10 — 0.45 as shown below. This value of A is 
larger than in the MROl-like model because the adopted mass ranges both for primary and 
secondary stars are narrower. 

ii) Mass ratio — In Figure [ini we show several functions of the mass ratio f{q) or the 



mass fraction /(//) (Note that q = 1712/ mi 



^ and ^ 



(ilq-^i = (1 — /^)^)- The function 
that is introduced in MROl is weighted for massive secondaries (i.e., the positive value of 
7, long-dashed line). Among the binary populati on synthesis, several fu nctions are adopted: 
/(g) oc (1 + q)^'^, i.e., constant /(/i) f dott ed line. iPortegies Zwartlll995l ). and constant /(<?), 
i.e., /(/i) oc 1/(1 — /i)^ (short-dashed line, Han et al. 19951 ). Observations of binary systems, 
however, suggest that the function is weighted for less-massive seconda ries, as in our function 
(solid line). The histog ram shows the observation of G-dwarfs by iDuauennov &: Mavor 
(1991). For B-type stars ( Shatsky fc Tokovinin 2002; Kouwenhoven et alT2005 ). the power- 
law distributions with the slopes of —0.33 to —0.5 are reported, and the random paring of 
stars is not supported. Note that the shape of /(g) is affected by selection bias (IHogeveen 
I992I ). and may depend on the density a nd composition of the stellar aggregate where the 
binaries have been formed (jKroupalll995l ). Our function is between these observations; the 
secondary IMF 0d is weighted for less-massive secondaries, and the slope is shallower than 
that of the primary IMF 0. 



In Figure [201 we show that the mass ratio function does not change our results. In 
other words, it is possible to obtain almost identical results with different sets of A and 
/(/i). For our SN la model, although MROl adopted 7 = —0.35, 7 = —1.35 is correct 
because the number function should be the mass function 0d multiplied by 1/m. Since 
/(/x) with 7 = —1.35 is weighted for low-mass secondaries, the relative contribution of the 
MS-I-WD systems is small, and A should be as large as 0.45. With the constant /(/i) or 
constant /(g), the contribution of the MS-I-WD systems is as large as that of the RG-I-WD 
systems, and thus A becomes 0.30 and 0.25, respectively. With 7 = 2 as in MROl, A = 0.20 
gives the best fit. 



iii) Binary IMF — The IMF is applied to the total mass of the binary in Eq.[T], while 
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it is applied to the mass of the primary star in Eq.[2]. This may correspond to different 
approaches in the formation of binaries: a common or individual formation. This treatment 
also does not change our results. If we adopt not 0(mB) but 0(mi), the normalization is 
changed by a factor. For 7 = —1.35, constant /(yu), constant f{q), and •y = 2, A becomes 
0.25, 0.15, 0.12, and 0.10, respectively. 

Needless to say, for the models that give the same [0/Fe]-[Fe/H] relation, the SN la 
rate history is almost the same. In the lower panels of Figures ITS] and [2U| the time evolution 
of SN la rate is shown. Only 0.05 dex difference is seen in the present SN la rate among 
the models with different binary mass fractions. The present SN la rate with 7 = —1.35 is 
0.1 dex larger than with 7 = 2. It seems not possible to put constraint on the binary mass 
fraction and binary IMF from the observed SN la rate because of the large errorbar. 

As a summary, once a suitable 6 or A is adopted, almost identical results can be produced 
with either of these two formulations. The most important factor is the mass ranges of 
companion stars, which changes the distribution function of SN la lifetimes. The total 
number of SNe la, b or A, are determined from the observational constraints in this paper, 
but niay be larger than what can be obtained with binary population synthesis models (e.g.. 



Maozll2008l ). Other assumptions, the mass ratio function and IMF, are not very important 
since these two effects can be cancelled out by the choice of the total number of SNe la. 
Observations of binary systems and the theory of binary formation could give a clue to solve 
this degeneracy, and may put more constraints on SN la models. 

Table 6. Test results and predictions of the SN la models. 



metallicity effect 


SD 

yes no 


MROl 


DD 


Lifetime [Gyr] 


0.1 - 20 peak at 0.1, 1 


~ 0.3 


~ 0.1 


[a/Fe]-[Fe/H] relation 





X 


X 


Present SN la rate in Es 








X 


Cosmic SN la rate history 


peak at 2 ~ 1 constant 


increase 


increase 


SNe la observed at z > 2 in Ss? 


few yes 


yes 


yes 
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6. Conclusions 

We construct a new SN la model, based on the SD scenario for the SN la progenitors, 
taking account of the metallicity dependences of the WD wind (K98) and the mass-stripping 
effect on the binary companion star (HKN08). The lifetime distribution of SNe la is deter- 
mined from the main-sequence mass range of the companion stars in the WD binary systems, 
and is given as a function of metallicity (Fig. [3]). The SN la rate in the systems with [Fe/H] 
< — 1, e.g., r > lOre in intermediate-mass galaxies, t > 50 kpc in dwarf galaxies, and high-z 
spiral galaxies, is supposed to be very small ( §2.3p . 

Our model naturally predicts that the SN la lifetime distribution spans a range of 
tia ~ 0.1 — 20 Gyr with the double peaks at ~ 0.1 and 1 Gyr reflecting the two systems of 
the companion stars: the MS+WD and RG+WD systems, respectively. With the chemical 
enrichment history in the Galactic disk, the fraction of the young population of SNe la is 
10% and 50% for < 0.1 and 1 Gyr, respectively (Fig. [TU]). In the galaxies with a initial 
star burst, the SN la rate shows a strong peak at t ~ 1.5 Gyr, which may generate the SN 
la driven galactic winds. 

We make a comparison of the lifetime distribution functions derived from different SN la 
models, and evaluate the results from the observational constraints as summarized in Table 
El 

(1) Contrary to the argument in Matteucci & Recchi (2001, MROl), our SN la models give 
better reproduction of the [a/Fe]-[Fe/H] relation in the solar neighborhood. In the DD and 
MROl-like models, the typical lifetimes of SNe la are ~ 0.1 and 0.3 Gyr, respectively, which 
results in the too early decrease in [a/Fe] at [Fe/H] ~ —2. If we do not include the metallicity 
effect in our model, [a/Fe] decreases too early because of the shortest lifetime, ~ 0.1 Gyr, of 
the MS+WD systems. In other words, the metallicity effect is more strongly required in the 
presence of the young population of SNe la, to be consistent with the chemical evolution of 
the solar neighborhood. 

(2) At the same [Fe/H], [Mn/Fe] starts to increase toward higher metallicity, because SNe 
la produce [Mn/Fe] > 0. This evolutionary change is started from the same [Fe/H] as the 
decreasing trends of [a/Fe]. [Zn/Fe] evolves differently for the different SNe la models, 
because of the combination of the different lifetime of SNe la and the metallicity effect on 
Zn yields of SNe 11. Further observations of these elements at [Fe/H] > — 1.5 is critically 
important in order to identify the SN la progenitors and to clarify the metallicity dependence 
on the nucleosynthesis yields. The lack of the similar [Mn/Fe] increase in dSphs suggests 
that those are not enriched by SNe la. 

(3) To explain the observed SN la in ellipticals at 2; = 0, the lifetimes of a sufficient fraction 
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of SNe la should be longer than ~ 10 Gyr, which is satisfied in our SD model and the MROl- 
like model, but not in the DD model. Our SN la models are also su ccessful in reproducin g 
the SN II, Ibc, and la rates with their dependence on the galaxy type (IMannucci et al.ll2006l ): 
This cannot be reproduced with the DD model. The large SN la rate in radio galaxies could 
be explained with the young population of the MS+WD systems in our model. 

(4) We also provide the cosmic supernova rate history as a composite of those in spirals 
and ellipticals. The cosmic SN la rate in comoving density shows a peak at 2 ~ 1 and 
decreases toward higher redshift in our model. In contrast, this rate increases toward higher 
redshifts in the DD and MROl-like models. Because of the metallicity effect, i.e., because 
of the lack of winds from WDs in the binary systems at [Fe/H] < —1.1, the SN la rate in 
the low- metallicity systems, e.g., high-z spiral galaxies, is supposed to be very small in our 
model. This metallicity effect will appear only around [Fe/H] ~ —1 (Fig. In contrast, 
the SN la rate in such systems is as high as present in the DD and MROl-like models. In our 
models, at 2; > 1, SNe la will be observed only in the systems that have evolved with a short 
timescale of chemical enrichment. This suggests that the evolution effect in the supernova 
cosmology can be small. 
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Fig. 1. — The distribution functions of the progenitor mass (upper panel) and hfetime (lower 
panel), i.e., the SN la rates of stars with the same ages, for K98 model with Z = 0.002 
(red solid line) and Z = .02 (red short-dashed line), MROl model (blue long-dashed line, 



Matteucci fc Recchil l200ll ). and DD model (green dot-dashed line, iTutukov fc Yungelson 
19941 ). 
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Fig. 2. — The distribution functions of the progenitor mass (upper panel) and hfetime 
(lower panel) for the metalhcity Z = 0.002 (blue), 0.004 (green), 0.02 (orange), and 0.05 
(red). The solid and dashed lines indicate the MS+WD and RG+WD systems, respectively, 
that correspond to the young and old population of SNe la. 
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Fig. 3. — The distribution functions of the progenitor mass (upper panel) and hfetime 
(lower panel), i.e., the SN la rates of stars with the same ages, for our new models with 
the metallicity Z = 0.002 (blue short-dashed line), Z = 0.004 (green solid line), Z = 0.02 
(orange dotted line), and Z = 0.05 (red long-dashed line). 
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Fig. 4. — The chemical evolution of the solar neighborhood: (a) the star formation rate, (b) 
the age-metallicity relation, and (c) the metallicity distribution function for different SN la 
models: our SD model (red solid line), SD model without metallicity effect (cyan dotted line), 
DD model (green short-dashed line), MROl model (blue lon g-dashed line), an d K98 model 
(magenta dot-dashed line). Observ ational data sources are: Matteuccil ( 1997), errorbar i 



m 



panel (a); lEdvardsson et al.l (Il993l ). filled circles in panels (b) and (c); IWyse fc Gilmore 



(Il995l ). open circles in panel (c) 
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Fig. 5. — The evolution of the elemental abundance ratios, [X/Fe]-[Fe/H] relations, in the 
solar neighborhood for different SN la models: our SD model (red solid line), DD model 
(green short-dashed line), and MROl model ( blue long-dashed line). Observational data 
sources are: For dis k stars, lEdvardsson et al.l (119931 ) . small open cir cles; thin disk s tars i n 
Bensby et al. J2003I . 2004), small filled circle s; accre tion component in lGratton et al.l (^2003 ). 
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Fig. 6. — Same as Figure [5] but for our SD model with (solid line) and without (dotted line) 
metallicity effect. 
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Fig. 7. — Same as Figure [5] but with different c; c = (0,1,3,5) (solid line, main model), 
c = (0,2,5,10) (short dashed line), and c = (2,3,5,10) (long-dashed line) respectively for 
Z = 0.002,0.004,0.02,0.05. 
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Fig. 8. — Same as Figure H] but for the SN la rate history in the solar neighborhood with 
different SN la models: our SD model (red solid line), SD model without metallicity effect 
(cyan dotted line), DD model (green short-dashed line), MROl model (blue long-dashed 
line), and K98 model (magenta dot-dashed line). The dot-dashe d line shows the SN II r ate. 
The dot shows the observed SN la rate in SOa-Sb type galaxies (jCappellaro et al. 1999 ). 
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Fig. 9. — The supernova rates against the metalhcity indicators, (a) the iron abundance, 
(b) the gas-phase metalhcity, and (c) the mean stellar metalhcity, with the star formation 
history in the solar neighborhood. 
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Fig. 10. — The cumulative function of SN la lifetimes in the solar neighborhood with our SD 
model at the galactic ages of t = 3 (clue dotted line), 5 (cyan short-dashed line), 10 (green 
long-dashed line), and 13 Gyr (red solid line). 
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Fig. 11. — The SN la rate history in elhpticals with different SN la models: our SD model 
(red solid line), DD model (green short-dashed line), and MROl model (blue long-dashed 
line). The dot-dash ed line shows the SN II rate. The dot shows the observed SN la rate in 
El-SO type galaxies (ICappellaro et al.lll999l ). 
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Fig. 12. — The time evolution of (a) the SFR, (b) the SN II and la rates, (c) the gas fraction 
per luminosity, and (d) the B-V color for four types of spirals (red dotted line, SOa/Sa; yellow 
solid line, Sab/Sb; green long-dashed line, Sbc/Sc; blue short-dashed line, Scd/Sd). In the 
panel (b), the upper and lower four lines are for the SN II and la rates per B-band luminosity 
(the s upernova units, SNu), respectively. Observational data sources are: ICappellaro et al. 
(|1999l). open circles (SN II) and filled circles (SN la) of SOa-Sb and Sbc-Sd in panel (b); 
Roberts fc Hayned (119941 ). crosses in panels (c) and (d). 
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Fig. 13. — The present supernova rates per K-band luminosity (upper panel) and per mass 
(lower panel) against the morphological tj^e of galaxies for SNe II (blue das hed line), la (red 



solid line), and Ibc (green dotted line). The observational data is taken from lMannucci et al. 
(120051 ) for SNe II (open circles), la (filled circles), and Ibc (crosses). 
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Fig. 14. — Same as Figure [131 but for different SN la models: our SD model with (red solid 
line) and without (cyan dotted line) metallicity effect, DD model (green short-dashed line), 
and MROl model (blue long-dashed line). 
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Fig. 17. — The stellar lifetime as a function of initial mass and met alii city. The long and short 
dashed-lines are for the analytic functions adopted in MROl and iHachisu. Kato &: Nomoto 
([20081), respectiyely. 




Fig. 18. — (a) The [0/Fe]-[Fe/H] relation in the solar neighborhood for different SN la mod- 
els. The solid line is for our original K98 model. The dotted line is calculated with MROl 
formula Eq. [1] with the appropreate parameters for K98 model (7 = —1.35,74 = 0.45). The 
dashed and long-dashed lines correspond to Model 1 and 3 in M ROl, respectively. Obser- 



vational data sources are: 



^or d isk stars, lEdvardsson et al.l (119931 ) . small op en circles; thin 



disk s tars in iBensby et al.l (120041). small filled circles ; accretion component in iGratton et al. 



(120031 ). triangles. For halo stars. ICayrel et al.l (j2004j ). large filled circles, (h ) The time evolu 



tion of SN la rate in SNu. The dot is for the observational data for SOa-Sb (iCappellaro et al. 
1999h . 
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Fig. 19. — The distribution function of mass ratios: (g) oc m~^^~^^^ with x = 0.35 (sohd 
hne, this work), constant f{fi) (dotted hne), constant /(g) (short-dashed hne), and /(/i) oc 
/ i'^ with 7 = 2 (lonR-dashe d hne, MROl). The histogram shows the observational data 
( jPuquennoy fc Mayoiill99ll ). 
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Fig. 20. — The same as Figure fTSl but for different SN la models. The sohd hue is calculated 
with MROl formula Eq. [1] with the appropreate parameters for K98 model. The dotted and 
short-dashed lines are for the different functions of mass ratios, constant /(/i) and constant 
/(g), respectively. The long-dashed line is for the same /(/i) as in MROl, but with our 
parameters of binary companions. 
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Table 1. The SN la rate of simple stellar population in [10 ^ Mq^] as a function of 

companion's mass rrid [Mq]. 
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6.0 























0.228 








0.142 


0.003 


6.5 
































0.088 
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Table 1 — Continued 











this work 








K9S 




MROl 


DD 




Z = 


0.002 


z = 


0.004 


z = 


0.02 


z = 


0.05 


[Fe/H] > 


-1.1 








^Ia,m 


^Ia,m 


^Ia,m 


^Ia,m 




"^la.m 






"^la,!!! 






^Ia,m 




RG 


MS 


RG 


MS 


RG 


MS 


RG 


MS 


RG 


MS 






7.0 
































0.049 





7.5 
































0.020 





8.0 
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Table 2. The SN la rate of simple stellar population in [10 ^ Mq^] as a function of time 

tu [Gyr]. 



this work 



Z = 0.002 Z = 0.004 Z = 0.02 Z = 0.05 





'Ma,t 


'Ma,t 


''la 


'Ma,t 


'Ma,t 


"-la 


'Ma.t 


'^Ia,t 




'Ma,t 


'^Ia,t 






MS 




J- LVJ 


MS 




RG 

J- LVJT 


MS 




RG 


MS 


45.15 








57.50 








84.76 








78.29 








27.16 








32.39 








48.76 








44.78 








17.08 








19.56 








29.17 








27.04 


0.021 





11.33 


0.048 





12.53 


0.039 





18.59 


0.026 





17.24 


0.086 





7.974 


0.219 





8.413 


0.181 





12.36 


0.119 





11.00 


0.124 





5.891 


0.305 





5.929 


0.253 





8.575 


0.166 





7.596 


0.193 





4.533 


0.401 





4.265 


0.338 





6.026 


0.229 





5.563 


0.271 





3.563 


0.487 





3.165 


0.480 





4.501 


0.327 





4.151 


0.374 





2.812 


0.587 





2.519 


0.649 





3.462 


0.441 





3.322 


0.491 





2.271 


0.558 





1.998 


0.776 





2.783 


0.562 





2.607 


0.557 





1.807 








1.629 


1.217 





2.234 


0.675 





2.082 


0.777 





1.502 








1.478 


2.083 





1.845 


0.840 





1.793 


1.109 





1.262 





0.099 


1.349 


2.127 


0.172 


1.540 


1.253 


0.101 


1.557 


1.418 


0.115 


1.063 





1.812 


1.225 


2.164 


2.801 


1.415 


2.129 


2.756 


1.412 


1.864 


2.413 


0.903 





2.678 


1.117 


1.758 


2.470 


1.305 


2.166 


2.803 


1.287 


1.597 


2.068 


0.836 





4.513 


0.979 





2.238 


1.198 


2.196 


2.842 


1.119 


1.421 


1.840 


0.777 





4.581 


0.864 





2.498 


1.105 


1.896 


2.455 


0.979 


1.599 


2.070 


0.720 





4.637 


0.766 





2.774 


0.981 


1.691 


2.189 


0.862 


1.790 


2.317 


0.669 





3.804 


0.683 





3.066 


0.876 


1.840 


2.382 


0.763 


1.973 


2.554 


0.596 





3.279 


0.612 





3.375 


0.783 


1.993 


2.580 


0.676 


2.164 


2.802 


0.533 





3.374 


0.550 





3.701 


0.702 


2.160 


2.796 


0.603 


2.372 


3.070 


0.475 








0.497 





4.045 


0.630 


2.330 


3.017 


0.538 


2.587 


3.349 


0.425 








0.450 





4.408 


0.568 


2.554 


3.306 


0.482 


2.852 


3.692 


0.386 








0.409 





4.789 


0.514 


2.789 


3.610 


0.434 


3.132 


4.055 


0.351 








0.374 





6.163 


0.467 


0.239 


4.431 


0.392 


3.912 


5.063 


0.268 








0.296 





6.864 


0.356 





4.899 


0.295 


4.371 


5.658 


0.212 








0.239 





8.623 


0.278 





6.398 


0.228 





7.479 


0.161 








0.184 





1.847 


0.206 





7.788 


0.167 





9.124 


0.135 








0.157 








0.171 





10.41 


0.137 





12.17 


0.114 








0.135 








0.144 





12.54 


0.114 





14.86 


0.094 








0.112 








0.116 





15.43 


0.090 





18.65 


0.075 








0.092 








0.093 





17.97 


0.071 





23.34 


0.062 








0.077 








0.075 








0.057 





9.631 


0.052 








0.065 








0.063 








0.048 








0.045 








0.056 








0.053 








0.040 









-54- 



Table 2 — Continued 



this work 

Z = 0.002 Z = 0.004 Z = 0.02 Z = 0.05 

^la T^Ia.,t T^Ia.,t ha. '^Ia,t "^la.t ^la T^Ia.,t T^Ia.,t ha. "^la.t '^Ia,t 

RG MS RG MS RG MS RG MS 



0.039 0.050 0.046 0.035 
0.034 0.044 0.041 0.030 
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Table 3. The SN la rate of simple stellar population in [10 ^ Mq^] as a function of time 

tu [Gyr]. 







K98 






MROl 


DD 


z 


= 0.002 




Z = 0.02 










T^ja, t 






^Ia,t 




^Ia,t 






RG 


MS 




RG 


MS 






45.15 








84.76 








0.003 





27.16 








48.76 








0.007 





17.08 








29.17 








0.017 





11.33 


0.062 





18.59 


0.034 





0.039 


0.006 


7.974 


0.284 





12.36 


0.154 





0.080 


0.017 


5.891 


0.395 





8.575 


0.215 





0.156 


0.039 


4.533 


0.519 





6.026 


0.297 





0.284 


0.087 


3.563 


0.630 





4.501 


0.423 





0.535 


0.168 


2.812 


0.760 





3.462 


0.571 





0.913 


0.302 


2.271 


0.722 





2.783 


0.592 





1.463 


0.449 


1.807 








2.234 








1.668 


0.624 


1.502 








1.845 








1.945 


0.832 


1.262 





0.345 


1.540 





0.351 


2.727 


1.090 


1.063 





6.289 


1.415 





9.565 


4.363 


1.239 


0.903 





9.295 


1.305 





9.730 


4.189 


1.399 


0.836 





15.66 


1.198 





9.865 


4.017 


1.589 


0.777 





15.90 


1.105 





8.520 


3.287 


1.794 


0.720 





16.09 


0.981 





7.596 


2.780 


2.144 


0.669 





13.20 


0.876 





8.267 


2.875 


2.543 


0.596 





11.38 


0.783 





8.956 


2.963 


3.012 


0.533 





11.71 


0.702 





9.293 


3.057 


3.543 


0.475 








0.630 








3.143 


4.165 


0.425 








0.568 








3.285 


4.868 


0.386 








0.514 








3.425 


5.657 


0.351 








0.467 








4.014 


6.542 


0.268 








0.356 








3.873 


9.824 


0.212 








0.278 








4.426 


14.24 


0.161 








0.206 








4.510 


22.33 


0.135 








0.171 








5.267 


25.27 


0.114 








0.144 








5.522 


27.56 


0.094 








0.116 








5.598 


12.86 


0.075 








0.093 








5.274 


1.004 


0.062 








0.075 








4.807 


0.086 


0.052 








0.063 








4.066 





0.045 








0.053 








2.949 
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Table 3 — Continued 







K98 






MROl 


DD 


z 


= 0.002 




Z 


= 0.02 








ha 






ha 


T^Ia,t 


T^Ia,t 


'^Ia,t 






RG 


MS 




RG 


MS 






0.039 








0.046 








1.532 





0.034 








0.041 















